Nonadiabatic scattering of a quantum particle in an inhomogenous magnetic field 



F. Dahlem\ F. Evers^, A. D. Mirlin^^^, D. G. Polyakov^, P. Wolfle^-^ 
^ Ecole Doctorate de Physique, Universite Joseph Fourier - Grenoble 1, 38042 Grenoble, France 
^Institut fur Nanotechnologie, Forschungszentrum Karlsruhe, 76021 Karlsruhe, Germany 
Institut fiir Theorie der Kondensierten Materie, Universitdt Karlsruhe, 76128 Karlsruhe, Germany 

(February 1, 2008) 

We investigate the quantum effects, in particular the Landau-level quantization, in the scattering 
of a particle the nonadiabatic classical dynamics of which is governed by an adiabatic invariant. As a 
relevant example, we study the scattering of a drifting particle on a magnetic barrier in the quantum 
limit where the cyclotron energy is much larger than a broadening of the Landau levels induced by the 
nonadiabatic transitions. We find that, despite the level quantization, the exponential suppression 
exp{—2Trd/S) (barrier width d, orbital shift per cyclotron revolution 5) of the root-mean-square 
transverse displacement experienced by the particle after the scattering is the same in the quantum 
and the classical regime. 

PACS numbers; 3.65.Nk 72.10.-d 73.43. Cd 



A well-known example of an adiabatic invariant is the 
magnetic flux enclosed by a trajectory of a classical elec- 
tron drifting in two dimensions in a magnetic field per- 
pendicular to the drift plane (i?xi?-drift). Although it is 
not a strictly conserved quantity, it implies strong con- 
straints for the motion of the particle. In particular, in 
a smooth potential the the guiding center drifts along 
equipotential lines because a shift from the equipotential 
would change the cyclotron radius, which violates adia- 
baticity. This principle underlies, e.g., the semiclassical 
theory of the integer quantum Hall effect . It provides 
a mechanism for "adiabatic localization" of electrons |^] 
and leads to a non-monotonic dependence of the width 
of the cyclotron resonance on the magnetic field j|] . 

In a series of recent papers |Q] , it has been emphasized 
that adiabatic invariance is crucial also for the dynam- 
ics of composite fermions, which are quasiparticlcs rele- 
vant to the fractional Quantum Hall effect (FQHE) . The 
composite fermions drift along lines of constant flux in 
an effective random magnetic field if the filling fraction 
sufficiently differs from 1/2, so that the average effective 
magnetic field is strong enough. Then almost all trajec- 
tories are closed and it is only due to weak violation of 
adiabatic invariance that the conductivity is non-zero . 

Because of its significance for the understanding of the 
phenomena mentioned above, we have recently formu- 
lated and solved a model problem: the scattering of a 
classical particle drifting in a magnetic field with homoge- 
nous gradient towards a magnetic impurity |^] . We found 
that the displacement of the guiding center in a scatter- 
ing event is suppressed by an exponentially small factor 
exp(— 27rd/(5), where 6 is the shift of the guiding center in 
a cyclotron period along the drift trajectory and d is the 
impurity size. In Quantum Hall devices d is determined 
by the (spacer) distance between the doping layer and 
the plane of the two-dimensional electron gas. Moreover, 
we also found an unusual geometric resonance which cor- 
responds to "self-commensurability" . It is not related to 



the geometry of the scatterer but to that of the unper- 
turbed trajectory: whenever the cyclotron diameter 2Rc 
and S are commensurate the scattering is strong. 

It is an important question whether the strong adia- 
batic constraint on the possible classical scattering pro- 
cesses is lifted by quantum mechanics. If so, then the 
adiabatic-localization scenario outlined above and the 
corresponding interpretation of experiments could be 
questionable. In order to resolve these issues, we analyze 
in the present paper the scattering problem for a drift- 
ing particle in the quantum regime. Like in the classical 
case, we calculate the root-mean-square (r.m.s.) trans- 
verse shift Ap. In the quantum regime, the Landau-level 
quantization is strong, so that a separation in real space 
of two drift trajectories with the same energy and cor- 
responding to different Landau-levels by far exceeds the 
classically calculated Ap. Although one might expect 
otherwise, we find that the r.m.s. non-adiabatic shift is 
given by the same formula as in the classical case. 

An analogous conclusion has been drawn in earlier 
work for the case of motion in a random potential [^j. 
The difference between our consideration here and that in 
Ref. Q is twofold: i) our one-impurity problem has well 
defined in- and outgoing states; ii) we consider an inho- 
mogeneous magnetic field (rather than scalar potential). 
An important advantage of the one-impurity scattering 
problem is that it allows for a controlled calculation of 
the nonadiabatic shift. This yields not only the expo- 
nential suppression but also the prefactor. Thus, we find 
an additional non-trivial feature, the above mentioned 
"self-commensurability" oscillations. 

Let us point out further relations of this work to yet 
another field of interest. Recently, a number of ex- 
periments have been dealing with semiconductor hybrid 
structures which consist of a quasi one-dimensional con- 
ducting channel with a magnetic barrier that modulates 
the electronic current ■ These experiments have been 
performed in weak magnetic fields. At sufficiently strong 
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fields the device may be considered an experimental re- 
alization of our model problem. In principle, this allows 
for an experimental check of our predictions by measur- 
ing the transverse current or the Hall voltage. We come 
back to this issue at the end of the paper. 

The Hamiltonian for the free drift motion of the model 
problem formulated above reads 



Ho = (zV + Ao)V2^ 



(1) 



We choose Landau's gauge with Aq = By(l + i]y/2l,0'^, 

where I = y^h/mujc is the magnetic length, uJc ~ eB j mc 
is the cyclotron frequency, 77 is a dimensionless param- 
eter characterizing the magnetic field gradient and we 
adopted units such that e = c = 1. Due to translational 
invariance the exact eigenstates are of the form ip{x) — 
exp{ikx)xn{y) / S'Hd the problem is effectively one- 
dimensional. The functions Xn,k{y)i with n = 0, 1, 2 . . ., 
are eigenfunctions of a quantum particle in a quartic po- 
tential 
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After introducing new coordinates ^ 
obtains the effective Hamiltonian 



(2) 

(y ~ yk)/lk one 
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where h = + mk/l, yk = {y/T+TrJkl ~ 1)1 /ij and 

Vk = vih/l)^- The exact eigenfunctions Xn,k of hk are 
not known analytically and have been calculated numer- 
ically 1^. The simplest approximation is to ignore the 
anharmonic terms altogether. This procedure is justified 
if the relative change of the magnetic field within the cy- 
clotron radius = \J1n -I- IZ^ is small: rjkRkl'^.lk <C 1, 
which is what we assume here. In other words, we con- 
sider the case where the quartic potential forms a 
double well with two well separated parabolic minima. 
Our focus is on particle energies not too large such that 
tunneling between the two wells is strongly suppressed 
and the shape of each is approximately parabolic. In this 
situation the electronic trajectories resemble a cyclotron 
orbit drifting along lines of constant magnetic field. 

We are then facing a well-defined scattering problem 
with free scattering states 



Xn,k ~ MO = HniO exp (-eV2) /NnM- 



(3) 



Here, _ff„ denotes the Hermite polynomials and Nn^k = 
(-y7r2"n! 1^) " is the normalization factor. To leading 
order in 77 and 1/n the spectrum is given by enivk) = 
fuvkin + 1/2) with u)k = u;c(l -I- vyk/l). 

The magnetic scatterer has a vector potential SAx{x.) 
which gives rise to a perturbation to the Hamiltonian Hq 



where we have dropped the second-order term. 

Because of energy conservation any scattering process 
from fc to fc', i. c., from yf. to j/j,/ is accompanied by a si- 
multaneous change in the Landau level quantum number 
n. In the parabolic approximation, one has 



Vk 



yk' = {k-k')f - {2Trf/d) {n'-n), 



(4) 



where I - 



l + T]Y/l, Rc = ly/{n + n' + l), r] = 

iVk + yfe)/2 and S = TT{rj/l)Rc . The prob- 
ability w{nk; n') of a scattering event from an incident 
state with quantum numbers n, k to an outgoing state 
with n' can be calculated by means of the golden rule: 

■w{nk;n') — P{nk; n'kf)/Ink 

with the scattering probability per time P{nk]n'k') — 
(27r//i)(5(e„.fc — en'k')\M {nk; n' k')\'^ and the incident cur- 
rent In k — {fi'Lx)~^denk/dk. Within the parabolic ap- 
proximation one finds 



w{nk;n') = \AI{nk;n'k')\ 



/;/m2 / ^^nkd^n'k' 



dk dk' 



where k' is determined by (^). The drift velocity is 
denk/dk=SkL0k/2Tr, where Sk denotes the shift in x- 
direction of the corresponding classical orbit within the 
cyclotron period, Sk—Tr{r]k/lk)R'i in the parabolic ap- 
proximation. The corresponding matrix elements are 



M{nk-n'k') = (xnfcK 



Ax{y) k + k' 



2m 



)SAx{q,y)\x n' k' 



where q = k—k' and SA^ {q, y) is the Fourier transform of 
the scattering potential with respect to the x-coordinate. 

Before we proceed with evaluating M{nk, n'k') let us 
discuss different relevant regimes. If the scattering prob- 
ability w{nk, n') as calculated above from the golden rule 
turns out to be large, w{nk,n') 3> 1, we are in the qua- 
siclassical regime. In this limit of a strong barrier, when 
the Born approximation is not valid, the classical results 
derived earlier [|j should be applicable. For a weaker bar- 
rier, w{nk, n') <^ 1, the system is in the quantum regime, 
which is of central interest for us in the present paper. 

It should be noted that the Born approximation 
matches the classical description at w{nk, n') ~ 1 di- 
rectly only if the parameter bo/Srj <C 1, where b is the 
characteristic magnetic field produced by the impurity. 
If bo/Sr] <^ 1, one can neglect the curving of the quasi- 
classical trajectories by the impurity field. If bo/Srj ^ 1, 
the Born approximation fails with increasing b when 
w{nk, n') is still small, so that there appears an inter- 
mediate regime between the Born and classical limits. In 
this regime the scattering probability w{nk, n') can be 
calculated by means of an eikonal approximation. The 
crossover to classics occurs when the exact scattering 
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probability is of order unity. For our purpose, however, 
it is sufficient to consider the Hmit bo/Sr] <^ 1, which we 
assume in what follows. ^ 

We will further assume the condition 27rl /5 <^ Rc 
(which can be equivalently rewritten as tjVtt^ ^ 1 and 
which ensures that the spacing Sy between the available 
states is small compared to the cyclotron radius). In par- 
ticular, this condition implies that the Landau level index 
is large, n ^ 1. In the opposite case rjVri? <C 1 (which 
can be classified as an ultra-quantum limit) the states 
with n' ^ n overlap exponentially weakly, which leads to 
an additional suppression of the scattering cross-section, 
unrelated to adiabaticity. 

Now we return to the evaluation of the matrix element 
M{nk,n'k'). Since it can be calculated only approxi- 
mately, we first need to analyze the impact of gauge in- 
variance. In order to do so, we consider a particular 
example, namely a magnetic field barrier with a shape 
independent of y, SAx{q,y) = {y + yo)b{q)- Since the 
scattering amplitude must be independent of the choice 
of the origin yo we should have 

{Xnk\A.Ay) -{k + k')/2\xn'k') = (5) 

whenever the states Xnk,Xn'k' have the same energies. 
Indeed, this property of the exact eigenstates follows 
from the identity {xnk\hk — hk'\Xn'k') = 0, immediately. 
Apriori, we cannot expect the parabolic approximation 
(^ to fulfill this identity and in fact it does not. To 
further illustrate this point, we introduce the magnetic 
field b{q,p) corresponding to the scatterer by inverting 
b{,q,y) — —dySAx{y) in Fourier space 

M{nk;n'k')^ J g b{q,p)/{~tp) 

{Xnk\{AAy)-ik + k')/2)e^Py\Xn'k'). 

Due to Eq. (^) we can be sure that the pole p = 
is canceled if the exact matrix element is used. In or- 
der to avoid spurious terms resulting from any approxi- 
mation that could spoil this cancellation, we employ an 
exact rewriting of the formula and replace exp(ipy) 
exp{ipy) - 1. 

Even in the parabolic approximation Xnk ~ 4>nk the 
matrix element cannot be evaluated analytically, since 
the incoming and outgoing states experience slightly dif- 
ferent magnetic fields so that Ik and Ik' differ from each 
other. However, the difference may be ignored if the con- 
dition n, n' ^ \n' — n\ is met. This is precisely the case 
of interest to us, since n,n' ^ 1 and since the relevant 
values of rt' — n are ±1, as will be explained below. 

Therefore we simplify the matrix element by letting 
Zfc w 1^1 w I. Dropping further terms small in rj, 1 /n we 
obtain (n' > n) 

1 fdp 



M(nk, nk') = — 
fnn' ('*) 



2mAf 

n'-UQ 



27r 



b{q,p) e'P''U„,{{q + ip)l/2) 
i2\K'\) (6) 



with the normaHzation factor M = (2" n!/2"n'!)~^/^. 
The expression for the other case n < n' is found by 
interchanging n' , n and replacing q —q in fnn' ■ 

Let us discuss our result by considering several rele- 
vant examples. First, we consider the barrier SAx{q, y) — 
yb{q), again. We will see that higher order transitions 
|n' — n| > 2 are strongly suppressed if the width of the 
barrier d substantially exceeds S. Since it is precisely 
the adiabatic limit of smooth scatterers, d ^ (5, we are 
mostly interested in, we specialize to the case n' — n + 1 
right away. Also, we focus on displacements small as 
compared to the cyclotron radius, yk — yk' < Rc, which 
allows us to use the asymptotic expansion of the Laguerre 
polynomials. This yields for the transition probability 



w(jik; n + 1) 



1 



2Trl 



b{qi) 



^2^08 [qiRc 
qil 



where qi = 2n/5. For further discussion, we list 
two examples: relevant to the FQHE is the form 
6i(x) = h^d^ + d^)~^/'^ . Its one-dimensional Fourier- 
transform is bi{q,y = 0) = 6o(i\/27r|g|dexp(— l^jd) , valid 
at \q\d ^ 1, which makes the prefactor in the transition 
probability. 



wUnk; n + 1) - AW^^e-^'^^'^ cos' (q^R, - J) (8) 



independent of r/ (here Wq — bo/mLj). The factor 
exp(— 2(7ic?) represents the exponential suppression of the 
transition probability due to adiabaticity. For the pro- 
cesses with |n — n'l > 1 this suppression is much stronger, 
w{nk,n') oc exp(—2qid\n — n'\), so that they can be 
neglected, as stated above. For the magnetotransport 
experiments another magnetic field configuration 

is relevant, with an average value equal to zero and 
two "bumps" of opposite sign separated by a distance 
a, i.e. b{x) — (3{x — a/2) — (3{x + a/2). We choose 
(3{x) = bo\n{{x' + d')/{x'' + {d+D)')) |o|, which implies 



-2qid 



X sinh2(giD) cos^ (qiRc - ^) e^'^^. (9) 

In addition to the adiabatic suppression, we observe an 
oscillating behavior as the slope rj is varied, with a geo- 
metrical resonance condition a = jo5,jo = ±1, ±2, ... at 
which the leading order scattering process is forbidden. 

Next, we turn to another example where the scatterer 
is of finite extent also in ^/-direction. We evaluate the in- 
tegral for the case of 6i, as defined above, in two different 
limits, again focussing on n' = n -I- 1. When Rc exceeds 
all length scales, d,Y <^ Rc, we obtain 
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-2qi 



(10) 



for the transition probability. In the other regime, where 
Rc<,d and y = 0, we recover the ear her result Eq. (||). 

We are in a position now to compare the classical and 
quantum mechanical results for the quantity Ap which 
is defined as the r.m.s. magnitude of the shift perpen- 
dicular to drift velocity after the scattering took place, 
Ap = ((5y)2>i/2 ^ {2nfr5){T.n'{^^-^'?^nk,n'Y'^. For 
the case d and F = we have 



Apia = 25/2,1 Wo yRcdyS e'^'^'cos (qiRc-n/i) , 

(11) 

whereas when d,Y <^ Rc we find 

Apib = 2^/2^ Wo (d^/S) e-«i'^ cos (qji^) (12) 



under the conditions rj^/n ^ 1 <C rjvn^ and d/S 3> 1. 
These results arc identical with the classical findings re- 
ported in Ref. |^ after the latter have been averaged 
over the initial conditions. Note, however, that the pic- 
ture of the scattering is totally different in the classical 
and the quantum regimes. In the former case, it is the 
typical magnitude of the adiabatic shift which is expo- 
nentially small 6y cx exp{—2TTd/S). On the other hand, 
in the quantum case the scattering shift is set by the 

_2 

spacing between the unperturbed states, Sy = 2ttI /S, 
while the scattering probability is exponentially sup- 
pressed w cx exp{—4-iTd/S). So, the agreement found be- 
tween the results for the r.m.s. shift is far from trivial. It 
is also worth mentioning that it is exactly the r.m.s. shift 
that determines the macroscopic diffusion coefficient in a 
system with a finite density of scatterers . 

It is known that in the absence of a magnetic field the 
transport scattering rate has the same value in the classi- 
cal and quantum regimes, provided the correlation length 
of the randomness is large, kpd ^ 1 ||ll[]. We find now 
an analogous result for the case of transport in strong 
magnetic fields, when the dynamics is governed by the 
adiabatic invariance. 

Next, we briefly discuss scattering events in which the 

-2 

displacement yk —ye = qI exceeds the cyclotron radius, 

_2 

ql > Rc ("ultra-quantum limit"). In other words, we 



focus on very small gradients rjVn^ < 1. Once again we 
consider the magnetic barrier independent of y. Then 
the scattering probability is 



w(nk, n') — — ^ 
27rl 



<rJ_ Hi) 

R^ mU 



2 / -2 \ 2(n'+n) 



xe 



(13) 



provided that \n' — n\ n,n' . In the typical case, where 
the magnetic scattering field has a Fourier transform 
b{q) oc e"^"^, the Gaussian factor exp(— g^J ) reflecting 



the shape of the eigenfunctions starts to dominate over 

_2 

the adiabatic factor exTp{—2qd) at ql w max(2d, i?c)- 

Up to now we have considered the motion of a particle 
drifting in a magnetic field with homogenous gradient. 
The same analysis can be done also for particles that 
perform an i?xi?-drift. The transition probability is now 
w{nk;n') ~ \M{nk;n'k')\'^ /v"^ where we have chosen the 
_E— field E = {Ex,0) and introduced the corresponding 
drift velocity = E^/mujc. The matrix element is given 
by Eq. (^. With the substitutions S = 27rvd /ujc and 
rj = the results (0)-(^^ derived thereof remain valid. 

We have mentioned already that in general the scat- 
tering to the left and to the right have different probabil- 
ities. In the case of the E'xB-drift the transverse current 
resulting from processes with n' = n ± 1 can be easily 
evaluated: 

ly = ricUJc {b{q) /muJcY cos{2qRc) /tt, 

where Ug denotes the particle density. Similar to the case 
of Ap the current becomes very small whenever the con- 
dition d:s> S is satisfied. Therefore we predict a dramatic 
decrease of the transverse current upon crossing over into 
the regime of adiabatic dynamics at the barrier. 

Finally, let us mention that the calculation above im- 
plies a constant density of states. However, for electrons 
in smooth disorder the density of states splits into sep- 
arate Landau bands for sufficiently large B. It remains 
to be seen what will be the implications of adiabaticity if 
the crossover to drift occurs at still larger B. We expect 
a pronounced signature of the crossover to adiabaticity 
in the amplitude and shape of the magnetoresistance os- 
cillations. 
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